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Quantum graphs are commonly used as models of complex quantum systems, for 
example molecules, networks of wires, and states of condensed matter. We con- 
sider quantum statistics for indistinguishable spinless particles on a graph, con- 
centrating on the simplest case of abelian statistics for two particles. In spite of 
the fact that graphs are locally one-dimensional, anyon statistics emerge in a gen- 
eralized form. A given graph may support a family of independent anyon phases 
associated with topologically inequivalcnt exchange processes. In addition, for suf- 
ficiently complex graphs, there appear new discrete-valued phases. Our analysis is 
simplified by considering combinatorial rather than metric graphs - equivalcntly, a 
many-particle tight-binding model. The results demonstrate that graphs provide an 
arena in which to study new manifestations of quantum statistics. Possible applica- 
tions include topological quantum computing, topological insulators, the fractional 
quantum Hall effect, superconductivity and molecular physics. 
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1. Introduction 

The quantum mechanical properties of a many-particle system depend profoundly 
on whether the particles are distinguishable or not, even if, in the former case, the 
distinguishing properties have little effect on the underlying classical mechanics. 
For this reason, the description of indistinguishable particles in quantum mechan- 
ics - quantum statistics - continues to be revisited, both to gain insight into its 
foundations (Duck & Sudarshan 1997, Berry 2008) as well as to predict and un- 
derstand new phenomena, including the fractional quantum Hall effect (Wilczek 
1990, Jain 2007), superconductivity (Wilczek 1990), topological quantum comput- 
ing (Nayak et al. 2008) and particle-like states of nonlinear field theories, where 
standard quantization procedures are not straightforward to apply (Finkelstein & 
Rubenstein 1968, Manton 2008). 

The standard treatment of quantum statistics takes as its starting point the 
quantum mechanical description of a single particle. The n-particle Hilbert space, 
H n , is then taken to be the tensor product of n copies of the one- particle Hilbert 
space, H. The fact that the particles are indistinguishable means that physically 
observable operators commute with permutations of the particle labels. This im- 
plies that Hn may be decomposed into physically distinct components characterized 
by different permutation symmetries. The symmetrization postulate restricts the 
physically realizable spaces to either the symmetric component, which obeys Bosc 
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statistics, or the antisymmetric component, which obeys Fermi statistics; compo- 
nents characterized by more complicated behaviour under permutations - paras- 
tatistics - are thus excluded. Finally, the spin-statistics relation determines which 
of the two allowed alternatives applies - Bosc statistics for integer spin, and Fermi 
statistics for half-odd-integer spin. 

There is another approach to quantum statistics which is based on the topology 
of the classical configuration space. It takes as its starting point the configuration 
space, M, of a single particle. The n-particle configuration space, C n (M), is then 
taken to be the cartesian product of n copies of the one-particle configuration space, 
with, crucially, coincident configurations excluded (no two particles can be in the 
same place) and permuted configurations identified. So defined, C n (M) may have 
nontrivial topology. In particular, curves along which the particles are exchanged 
(some particles returning to where they started, others ending up where another 
started) are regarded as closed curves (as permuted configurations are identified), 
and these exchange curves cannot be continuously contracted to a single point. If 
such curves exist, then C n (M) is multiply connected. 

Quantum mechanics on a multiply connected configuration space C, whether C 
describes n particles or just one, allows for an additional freedom in the quantum 
description, namely the choice of a representation of the fundamental group, tti(C), 
of C . The (single-particle) Aharonov-Bohm effect provides a familiar example: the 
configuration space C for a particle with charge q in the presence of an impenetra- 
ble solenoid, say along the z-axis, may be taken to be R 3 with the solenoid - an 
infinite cylinder - removed. C is multiply connected; closed curves may be classi- 
fied according to the number of times they wind around the excluded cylinder. To 
the canonical momentum operator, p = — iV, we can add a vector potential, A, 
whose curl vanishes away from the cylinder, so that the commutation relations are 
preserved, as are the classical equations of motion - there are no physically acces- 
sible magnetic fields. However, the line integral of A around the cylinder need not 
vanish, and may be realized by a magnetic flux <fr inside the cylinder, measured in 
units of the flux quantum qh/c. The flux determines a representation of 711(C) that 
assigns to a closed curve with m windings the phase factor exp(im^>). The value of 
the flux modulo 2tt - cquivalcntly, the choice of a representation of ni(C) - affects 
the quantum mechanical properties, although observable consequences vanish in 
the classical limit. 

More generally, on a multiply-connected configuration space C, the momentum 
is described by a covariant derivative which, in the absence of external fields, has 
zero curvature (so that the classical mechanics is unaffected) but which may en- 
gender nontrivial holonomy. The holonomy corresponds to a representation of the 
fundamental group m(C). The representation could be one-dimensional and there- 
fore abelian, assigning (commuting) phase factors to noncontractible closed curves, 
or higher-dimensional, assigning dxd unitary matrices (in general, non-commuting) 
to noncontractible closed curves. The unitary matrices act on the wavefunctions, 
which are taken to have al components (more generally, wavefunctions are taken to 
be sections of a <i-dimcnsional complex vector bundle) . 

Taking C to be the configuration space C n (M) of n indistinguishable parti- 
cles, Leinaas and Myrheim (1977) showed that representations of the fundamental 
group TTi(C n (M)) determine the possible realizations of quantum statistics. Simi- 
lar conclusions were reached by Laidlaw and de Witt (1971) using path integrals, 
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an approach that has been applied in a variety of settings, including cases where 
canonical quantization is difficult to carry out (Finkelstein and Rubenstein 1968, 
Balachandran et al. 1993, Manton 2008). 

If iri(C n (M)) coincides with the symmetric group, S n , then the topological ap- 
proach to quantum statistics yields the same possibilities as does the standard one, 
namely Bose or Fermi statistics, or, if one gives up the symmetrization postulate, 
parastatistics, which are associated with higher-dimensional representations of S n . 
This is the case for M = M. d with d > 2. As is well known, the situation is different 
for M = R 2 , i.e. for particles in the plane. The fundamental group 7Ti(C n (R 2 )) is 
the braid group B n (M. 2 ), which has infinitely many elements. The characterization 
of unitary representations of the braid group is not complete (see, e.g., Birman 
& Brendle 2004), but its one-dimensional (abelian) representations include a new 
type of statistics - anyon statistics (Wilczek 1990). Closed curves on which two of 
the particles wind around each other m/2 times (odd values of m correspond to an 
exchange of positions) are assigned a phase factor exp(imo;), independently of any 
phases associated with external gauge fields. The special values a = and a = it 
correspond to Bose and Fermi statistics respectively. 

Anyon statistics have been found to provide deep insight into phenomena involv- 
ing strongly interacting electrons, for example the fractional quantum Hall effect 
(Laughlin 1983). In the composite fermion theory (Jain 2007), for example, the 
collective excitations of the many-electron system are approximately described by 
an effective Hamiltonian in which they obey anyon statistics with phase a related 
to the fractional Hall conductance. 

The situation is different again for M = R, i.e. particles on the line. Here the 
fundamental group -k\ (C„(R 2 )) is trivial; there are no noncontractiblc closed curves, 
as it is not possible to exchange particles without one passing through another. 
Thus, the topological approach predicts only Bose statistics (but see Leinaas and 
Myrheim (1977) for yet another approach to quantum statistics in one dimension). 

One motivation for this paper is to study quantum statistics for particles on net- 
works of one-dimensional wires, or metric graphs. A metric graph T is a collection 
of nodes, or vertices, connected by edges, or one-dimensional intervals, of specified 
lengths. Single-particle quantum mechanics on metric graphs has been extensively 
studied. The model was originally introduced to describe free electrons on chemi- 
cal bonds (Pauling 1936). Some current applications include superconductivity (de 
Gennes 1981), quantum chaos (Kottos and Smilansky 1997, Keating 2008), and 
Anderson localization (Aizcnman et al. 2006). The particle is described by a set of 
wavefunctions, ^ t {x e ), one for each edge of T. Each wavefunction is acted on by a 
Hamiltonian, H e . For example, H e could be a Schrodinger operator with a gauge 
potential, i.e. of the form ^(—id/dx e — A e (x € )) 2 + V e (x e ). 

The boundary conditions at the vertices of the quantized graph is a subtle point 
in the theory. For free particles, i.e. H e = —^d 2 /dx 2 , a simple choice motivated by 
physical considerations is to require that i) wavefunctions are continuous across 
vertices and ii) the sum of the outgoing derivatives at each vertex should vanish 
(Neumann-like boundary conditions) . A complete characterization of boundary con- 
ditions that render the free-particle Hamiltonian self-adjoint has been given in terms 
of the Lagrangian planes of a certain complex symplectic boundary form (Kostrykin 
and Schrader 1999, Kuchmcnt 2004). For more general operators, the classification 
of self-adjoint boundary conditions remains an open problem (see Bolte and Har- 
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rison (2003) for a corresponding classification scheme for the Dirac operator). For 
a given graph model, one can in principle derive physically appropriate boundary 
conditions from the limiting behaviour of a more realistic two- or three-dimensional 
model, but the analysis can be quite involved. 

As a metric graph is locally one-dimensional, one might imagine that quantum 
statistics on a metric graph would be limited to Bose or Fermi statistics. This turns 
out not to be the case, as was previously demonstrated in particular examples by 
Balachandran & Ercolessi (1992) (discussed in Section BBf . However, a systematic 
treatment of quantum statistics for general metric graphs has not been given. One 
difficulty is that the quantization of many-particle metric graphs poses analytical 
challenges beyond those encountered in the one-particle theory. Besides boundary 
conditions for a single particle at a vertex, one has also to specify boundary condi- 
tions for coincidences of two or more particles on edges and at vertices that render 
the many-particle Hamiltonian self-adjoint. 

Here we circumvent these difficulties while preserving the essential topological 
features by considering the simpler problem of indistinguishable particles on a com- 
binatorial graph G. A combinatorial graph consists of a set of vertices together with 
a specification of pairs of vertices that are connected by edges. The edges themselves 
do not contain any points, so the configuration space for the particle consists of the 
vertices alone, and hence is zero- rather than one-dimensional. Quantum mechanics 
on a combinatorial graph is a tight-binding model on the vertices. Thus, quantum 
statistics on combinatorial graphs is an interesting problem in its own right. 

It turns out that the n-particle configuration space, C„(G), can also be regarded 
as a combinatorial graph, and hence, in contrast to a many-particle metric graph, 
is easily quantized. Following the topological approach, we can characterize the 
possible quantum statistics in terms of representations of the fundamental group of 
C n {G). As C n (G) is a discrete space, it is not immediately apparent what is meant 
by its fundamental group, but one can define a combinatorial version that coincides 
with its metric- graph counterpart. 

For simplicity, we consider here two particles on a combinatorial graph, and 
restrict our attention to abelian (one-dimensional) representations of the funda- 
mental group. It turns out that a given graph may support a family of independent 
anyon phases associated with topologically inequivalcnt exchange processes. In ad- 
dition, for sufficiently complex graphs, there appear new discrete-valued phases. 
(Mathematically, these discrete phases correspond to representations of the torsion 
component of the abelianizcd fundamental group, or, equivalcntly, the first homol- 
ogy group of C n (G).) The results extend straightforwardly to abelian statistics 
for n > 2 particles. The extension to nonabelian statistics will be discussed in a 
forthcoming paper. 

The paper is organized as follows. In Section [5] we present quantum mechanics 
on a combinatorial graph as a tight-binding model. We introduce gauge potentials, 
the analogue of vector potentials, which assign phase factors to the edges of the 
graph. These phase factors are then incorporated into the off-diagonal matrix el- 
ements (transition amplitudes) of the Hamiltonian matrix. Gauge potentials are 
determined, up to a choice of gauge, by the accumulated phases along, or fluxes 
through, the cycles of the graph. 

In Section [3] we consider two indistinguishable particles on a combinatorial 
graph. The two-particle configuration space can itself be regarded as a combinatorial 
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graph (usually larger than the original graph) , that can be quantized following the 
general prescription of Section [2J We characterize, up to a choice of gauge, topo- 
logical gauge potentials on the two-particle graph, which determine the quantum 
statistics. Topological gauge potentials have the property that the only cycles with 
nonzero flux correspond to nontrivial closed curves in the metric setting. They are 
parameterized by a set of free statistics phases, which range between and 2n, as 
well as, in some cases, a set of discrete statistics phases whose values are constrained 
to certain rational multiples of 2n. A certain subset of the free statistics phases may 
be attributed to Aharonov-Bohm flux lines threading the one-particle graph; the 
remainder describe many-body effects f Section I3E|) . Bose and Fermi statistics, as 
understood in the conventional approach, may be recovered from particular choices 
of topological gauge potential (Section Our results are illustrated by a number 
of examples in Section 21 and in the concluding remarks we consider perspectives 
for further investigations and possible applications. 

2. Quantum mechanics on combinatorial graphs 

(a) Combinatorial graphs 

A combinatorial graph G consists of a set V = {1, . . . , v} of sites, or vertices, 
labeled 1 through v, which may be connected by bonds, or edges. It is convenient 
to describe the connectivity of the graph, i.e. its edges, by a v x v adjacency matrix, 
A, whose (j, k) th entry gives the number of edges from vertex j to vertex k. We 
assume that G is undirected, i.e. Ajk = A^j . We write j ~ k to indicate that j and 
k are connected by an edge. The number of edges in the graph, e, is then given by 
e = 5 J2jk Ajfc- Edges will be labeled by Greek indices which take values between 
1 and e. Given an edge e, we let e< and e> denote its vertices of lower and higher 
index respectively. Sometimes it will be useful to assign an orientation to an edge; 
our convention will be that positive orientation corresponds to going from e< to e>. 

We also assume that G is simple, so that there is at most one edge between any 
two vertices, and that no vertex is connected to itself, i.e. Ajk is equal to zero or one, 
and Ajj is equal to zero. Any graph can be made simple by introducing additional 
vertices on its multiply-connecting and self-connecting edges. The valency, or degree, 
Vj, of the vertex j is the number of edges connected to j, i.e. Vj = J^k Ajk- Finally, 
we assume that G is connected, i.e. every pair of vertices is connected by some 
sequence of edges. In terms of the adjacency matrix, this means that, for any given 
vertices j and k, (A n )jk ^ for some n. 

While our focus is on combinatorial graphs, some aspects of our treatment may 
be motivated by considering metric graphs associated to a combinatorial graph. 
Given a combinatorial graph, G, we can associate to it a metric graph, T, by as- 
signing a length L t > to each edge e of G. On T, e is regarded as an interval 
[0, L e ], with identified with e< and L e identified with e>. One can define continu- 
ous curves on T, and a metric is obtained by taking the distance between two points 
to be the minimum length of continuous curves joining the points. One can consider 
closed curves based at a vertex * and define the fundamental group 7Ti(r, *) in the 
usual way. 

Next we recall some basic facts about the topology of combinatorial graphs 
(see, e.g., Hatcher 2001), which will play a role in their quantization. A path 
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V = (iojii) • ■ • i jn) on a combinatorial graph G is a sequence of vertices in which 
consecutive vertices are connected by edges, i.e. j r ~ jV+i (thus, consecutive ver- 
tices must be distinct). The length of a path is the number of edges along it. A 
single vertex, j, may be regarded as a path of zero length, and an edge, e, regarded 
as the path (e<, e>). Given a path p = (jo, ji, . . . , j n ), we define the inverse of p, de- 
noted p _1 , to be the path (j„,j n ~i, ■ • • , jo)- If p = (jo, ■ • • , jn) and g = (fc , ...,k r ) 
are two paths such that the last vertex of p coincides with the first vertex of q, 
we define the concatenation, or product, of p and q, denoted pq, to be the path 

(jo, • • ■ j jn, &1) • ■ • j k r ). 

The path qpp~ 1 r describes the product of paths q and r with an intervening 
retracing of the path p. We want to regard paths that differ by retracings of inter- 
mediate components as being the same. Thus, we introduce the equivalence relation 
qpp~ 1 r = qr. Paths that are equivalent to a zero-length path, i.e. to their initial 
vertex, are called self-retracing. 

A cycle is a path c = (jo, ji, . . . , j n ) that is closed, i.e. j n — jo. A cycle is 
primitive if its vertices, apart from the first and last, are all distinct. The set of 
cycles on a combinatorial graph can be characterized with the aid of a spanning tree, 
which we define next. A tree is a combinatorial graph whose only cycles are self- 
retracing. A subgraph of a combinatorial graph G is a combinatorial graph whose 
vertices and edges are subsets of the vertices and edges of G. A spanning tree, T, 
of G is a subgraph that is a connected tree containing all the vertices of G. Thus, 
given any two vertices of G, there is a path on T, unique up to retracings, that joins 
them. G has at least one spanning tree, and, unless G is itself a tree, more than one. 
(An iterative algorithm for constructing a spanning tree is to remove an edge from 
a primitive cycle of G until no primitive cycles remain). Clearly, any spanning tree 
of G has v — 1 edges. We let / denote the number of edges of G not in a spanning 
tree, so / is one minus the Euler characteristic of G, 

f = e-(v-l). (2.1) 

Let T denote a spanning tree of G and * a vertex of G. Let us label the edges of 
G that are not in T by an index <f>, 1 < <f> < f. Let c$(*) denote the cycle obtained 
by proceeding along the (unique) path with no retracings from * to <j) < on T , then 
from (f> < to 0> along qb, and finally from </>> back to * along the (unique) path with 
no retracings on T. We call c^(*) a fundamental cycle. The definition depends on 
the choice of T and *, but the dependence on * is easily accounted for, as 

c^(**) =pc (*)p-\ (2.2) 

where p is the (unique) path on T from ** to *. An arbitrary cycle c beginning and 
ending at * can be expressed, up to retracings, as a product of fundamental cycles 
and their inverses, i.e. 

<£(*)••.(&(*), (2.3) 

where sj — ±1. (More formally, the set of cycles based at * modulo retracings forms 
a group, the combinatorial fundamental group 7rf (G, *). 7rp(G, *) is a free group on 
/ generators ci(*), . . . , c/(*) and is isomorphic to the fundamental group, 7Ti(r, *), 
of a metric graph, T, associated to G.) 
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(6) Quantization 

We regard the set of vertices V of a combinatorial graph G as a configuration 
space. V might be the configuration space of a single particle, but we shall not 
restrict ourselves to this point of view. Indeed, in Section [3] we will consider com- 
binatorial graphs whose vertices represent configurations of two particles. We take 
quantum mechanics on G to be given by a tight-binding model on the set of ver- 
tices, in which short-time transitions arc allowed only between vertices connected 
by edges. The Hilbcrt space is <C" , with basis vectors \j). 1 < j < v, describing 
states in which the system is localized at one of the vertices. Dynamics is given by 
the Schrodingcr equation, i\ip) = H\ip), where the Hamiltonian H is a v x v her- 
mitian matrix with nonzero off-diagonal entries only between connected vertices, 
i.e. 

H jk =0i{j^k and j * k. (2.4) 

One example of a Hamiltonian is (minus) the combinatorial Laplacian, or the dis- 
crete kinetic energy, H = A — D, where D is given by Djk = VjSjk- 

(c) Gauge potentials 

A gauge potential Q on G is a v x v real antisymmetric matrix such that Qjk 
vanishes if j oo k. Given a path p = (jo, ■ ■ ■ ,j n ) on G, we define 

n-i 

Q(j>) = J2Q jrjr+1 . (2.5) 

r=0 

Clearly, if p and p' differ by retracings, then 

n(p) = Q{p'), (2.6) 

and if p and q can be concatenated, then 

ft{pq) = n(p) + n(q). (2.7) 

For a cycle c, we refer to 51(c) as the flux of fl through c. From (|2.2[) and (|2.6p . it 
follows that the flux through a fundamental cycle c^,(*) is independent of *; hence 
we write f2(c^), omitting * from the notation. For an arbitrary cycle c, it follows 
from (gZZl) that 

r 

0(c) = ^(-l)^0( % ), (2.8) 

3=1 

where <j>j and Sj are given by (|2.3[) . Thus, all fluxes are determined by the fluxes 
through a set of fundamental cycles. 

We incorporate a gauge potential Q into a Hamiltonian H by multiplying the 
transition amplitudes from j to k by the phase factors exp (iCljk). The new Hamil- 
tonian is given by 

H? k =H jkex p(iSl jk ) (2.9) 

Clearly, H n is hermitian and satisfies (|2.4|) (see Oren et al. 2009 for a similar 
construction) . 
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We can motivate the prescription (|2.9|) by making the following analogy with 
vector potentials on a metric graph. We note that the non-diagonal part of the 
Hamiltonian can be expressed as a real linear combination of the v(v— 1) transition 
Hamiltonians 

S {j ,k) = \j)(k\ + \k)(j\, A m =i(\j)(k\-\k)(j\), (2.10) 

where j ~ k, and, for definitcness, we take j < k. Let e denote the edge between 
j and k, and consider a metric graph where e corresponds to the interval [0, L e ]. 
Let Vt' e (x e ) denote the wavefunction along e on the quantized metric graph, and 
let us formally regard \j) and \k) as vertex states localized at x t = and x c = L f 
respectively, so that (i|* e ) = *e(0) and (k\^ e ) = ^ e (L e ). Then, formally, \k) = 
T e \j), where T e = e~xp(—ip e L e ) is the unitary translation by a distance L e along e 
generated by the momentum operator p e = —id/dx e . Thus we can rewrite (|2.10[) as 

%fc) = PjT} + , A [m = i (P,Tt - T e Pj) , (2.11) 

where Pj denotes the projector |j)01- Now suppose we introduce a vector potential 
on the edge e, replacing p e by p e — A e (x e ), and make the substitution 

T e -> cxp f-ij \p t - A e {x t ))alx\ = e l0 ^T e , (2.12) 



where Sljk is the integral of A e along e. Making this same substitution in (|2.11|) . we 

and ^8,«. 



obtain expressions for the transformed transition Hamiltonians, . » and Af 



that are equivalent to the prescription 

A gauge potential is trivial if it is of the form 

<U \' h «<--^>^ * 1 (2.13) 
I 0, otherwise, 

where 9 is a u-tuple of phases and M an (antisymmetric) integer matrix. For a 
trivial gauge potential, H n is given by H n = UHU\ where U is the diagonal 
unitary matrix Ujk = exp(i9j)8jk that generates the gauge transformation \j) n- 
exp(i9j)\j). The terminology is in analogy with quantum mechanics on M 3 , where 
a vector potential A(r) is trivial if it is a gradient V#, in which case it is induced 
by the gauge transformation 'f(r) i-> e ie( - r ' ) ^'(r). 

In analogy with the fact that a trivial vector potential has vanishing curl, we 
have the following characterization of trivial gauge potentials on G (we omit the 
argument, which is not difficult): 

n is trivial if and only if f2(c) is an integer multiple of 2n for every cycle c. 

(2.14) 

It follows from (|2.14l) and (|2.8|) that a gauge potential il is determined up to a gauge 
transformation by its fluxes through a set of fundamental cycles. Indeed, given a 
set of fluxes u>$, we can construct a gauge potential for which O(c^) = w^, as 
follows: if j and k are vertices of an edge in T, we take 

fl jk = fl kj = 0. (2.15) 

For an edge <f> not in T with vertices j — <j) < and k = 0> , we take 

Cljk = -flkj = w<f,. (2.16) 
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3. Two particles 

(a) Quantization 

Given a one-particle configuration space, X , there is a standard construction for 
the configuration space, C 2 (X), of two indistinguishable particles on X (the con- 
struction generalizes straightforwardly to more than two particles). C2(X) consists 
of unordered pairs of distinct points of X, i.e. 

C 2 (X) = {XxX-A 2 (X)}/S 2 , (3.1) 

where A 2 (X) = {(x,x)} denotes the coincident two-particle configurations, which 
are excluded, and S 2 denotes the symmetric group, whose single nontrivial element 
- exchange - acts on X x X according to (x, y) i— > (y, x). 

Given a combinatorial graph G, we shall now regard the set of its vertices, 
V = {1, ...,«}, as a one-particle configuration space. Hence, we take G 2 := C 2 (V) 
to be the configuration space for two indistinguishable particles on G. Depending 
on context, we denote configurations in G 2 either by unordered pairs of vertices 
{j, k}, or by ordered pairs (J, k) with j < k. 

We may regard G 2 as a new combinatorial graph. Nodes (or G2-vertices) (j, I) 
and (k,m) are taken to be connected by an edge if they have one G- vertex in 
common while their other G-vertices are connected by an edge of G. Equivalcntly, 
moving along an edge of G 2 corresponds to keeping one particle fixed at a vertex of 
G while moving the other along an edge of G. The adjacency matrix of G 2 , denoted 
A 2 , is given by 

^■2;{j,t),{k,m) = ^jk^-lm + 8jm.Alk + AjkSlm + Aj m Slk, (3.2) 

where A is the adjacency matrix of G. (To avoid awkward notation, we will not 
relabel the vertices of G 2 by a single index.) The number of vertices and edges of 
G 2 , denoted e 2 and v 2 respectively, are given by 

V2 = v(v-l)/2, e 2 = e(v-2). (3.3) 

We define quantum mechanics on G 2 following the prescription of Section [20 
The Hilbert space is C 2 with basis vectors \jl), j < I, representing states where 
one particle is at site j and the other at site I. Two-particle Hamiltonians, de- 
noted H 2 , are ^-dimensional hermitian matrices H 2 -ji^km, where j < I and k < to, 
whose off-diagonal elements, i.e., elements for which (J,l) ^ (k,m), vanish when- 
ever A 2: ( i j ! ;).(fe. m ) vanishes. Thus, short-time transitions generated by H 2 involve 
one particle making an allowed transition on G while the other particle remains 
fixed. 

Given a one-particle Hamiltonian H on G, we can construct a two-particle 
Hamiltonian on G 2 according to 

(jllHZlkm) = 8 jk H lm + aS jm Hi k + aSi k H jm + 5i m H jk , j <l,k <m, (3.4) 

where a = ±1 (one can check that satisfies (|2.4[) ). As discussed in Section l3fel 
below, a = — 1 describes nonintcracting fcrmions, while a = +1 describes hard- 
core (and therefore interacting) bosons which are prevented from occupying the 
same site. A simple illustration (two free particles on a linear graph) is given in 
Section PBfl 
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(b) Topological gauge potentials 

Gauge potentials can be introduced on a two-particle combinatorial graph fol- 
lowing the general prescription of Section Qfc] Quantum statistics is described by 
a subset of these, which we call topological gauge potentials. As discussed below, 
topological gauge potentials correspond to gauge potentials on a two-particle met- 
ric graph whose fluxes through all contractiblc closed curves vanish (modulo 2tt). 
Thus, their effects are purely quantum mechanical, vanishing in the classical limit. 

Let G be a combinatorial graph, and G2 the corresponding two-particle combi- 
natorial graph. Let e and <f> denote disjoint edges of G, i.e. e and cf> have no vertices 
in common. Let denote the cycle on G2 given by 



e,<t> = ({e<,0<},{>>, </><},{>>, </>>}, {><,</>>},{><, </><}) 



(3.5) 



That is, along c e ^, the particles move in alternating steps back and forth along e 
and 4>. We say that is metrically contractiblc The reason is that c e ^ corresponds 
to a loop 7^0 on C2(r), the two-particle configuration space for a metric graph T 
associated with G, which can be continuously contracted to a point - see Figure [TJ 
(These considerations may be made more formal. Let ** denote a vertex of the 



V 



Do 




a) 



b) 



c) 



Figure 1. a) The cycle c E ,0 and b) corresponding loop 75,^ in the two-particle 
configuration space of the metric graph Y. c) Contraction of 7 £[ 



two-particle configuration space. One can show that 7Ti(C2(r, **)) is isomorphic to 
the quotient of n c (G2,**) by the subgroup K(G2,**) generated by cycles of the 
form pc £i 0p _1 , where e and tfi are disjoint and p is a path from ** to {e<, 0<}.) 

Let £^2 denote a gauge potential on G2. We say that H2 is a topological gauge 
potential if its flux through every metrically contractible cycle vanishes modulo 27r, 
i.e. 

^(ce,^) = mod 27T, for e, <fi disjoint. (3-6) 

(More formally, topological gauge potentials may be identified with one-dimensional 
representations of ^(C^F), **).) 

Informally, gauge potentials that are not topological may be understood to 
introduce additional forces into the underlying classical dynamics, as the following 
heuristic argument, based on the analogy with metric graphs, demonstrates: when 
one particle is on the edge e and the other on 0, the free-particle Hamiltonian for the 
two- particle metric graph is given by \(j>\ If we introduce a gauge potential, 

replacing p e by p e — A e (x e ,y c j > ) and p^ by p^ — A^{x e , y^), the equations of motion 
become 

x e = B(x e ,y0)y , y = -B(x e , y<j,)x t , (3.7) 
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where the gauge field, B, is given by dA^/dx^ — dA t jdy^. If B ^ (analogous 
to ^(ce.c/,) mod 27r), the classical motion is no longer free. For combinatorial 
graphs, some numerical evidence for the effects of these gauge forces is given in 
Section 0|] 

(c) Characterization of topological gauge potentials 

Given a combinatorial graph G, we obtain below an explicit parametrization of 
the topological gauge potentials on the two-particle graph, Gi . The parametrization 
is given in terms of a set of fluxes, or phases. 

Let T2 denote a spanning tree of G2 ■ As shown in Section [^J given a gauge 
potential O2 on G2 (not necessarily topological), we can choose a gauge so that O2 
vanishes on the edges of T2. Let us label the edges of G2 that are not in T2 by an 
index (j>2, and let c^ 2 denote the corresponding fundamental cycles. The number of 
such cycles, which we denote by f 2 , is given by (cf. (|2.1|l ) 

h = e 2 - {v 2 - 1) = e(v - 2) - v(v - l)/2 + 1. (3.8) 

Let W0 2 denote the flux of O2 through c^ 2 . Then 0,2 is determined by u) — (<jJ\, . . . , w/ 2 ). 

The conditions 02(c e ,^) = mod 2ir on topological gauge potentials (cf. t\'3.Q]i ) 
can be expressed as linear relations on the components of uj. It will be convenient 
to label pairs of disjoint edges (e, 4>) by a single index a. The number of such pairs, 
which we denote by 172, is just the difference between the number of all pairs of 
edges of G and the number of pairs of edges that share a vertex, so that 

V 

.92= 2-e(e-l)-±5>>,-l). (3.9) 

Therefore, the condition for 0,2 to be a topological gauge potential can be written 

as 

i?-w = 27rn, (3.10) 

where R is a 32 X fi integer matrix and n an arbitrary ^-dimensional integer 
vector. The specific form of R depends on the choice of fundamental cycles. We 
note that the rows of R may be linearly dependent (this turns out to be the case 
if, for example, G contains two or more disjoint cycles). 

The system (|3.10[) can be solved by expressing R in Smith normal form (see, 
e.g., Dummit & Foote 2003). We write 

R = PDQ, (3.11) 

where P and Q are integer matrices with integer inverses of dimensions 52 and fa 
respectively, and D is a 52 X f% nonnegative diagonal integer matrix. The number 
of nonzero diagonal elements of D is given by the rank of R, denoted r, and D may 
be chosen so that its first r diagonal elements are nonzero with Djj divisible by 
Dj-i j—i for 1 < j < r. The nonzero diagonal elements of D are called the divisors 
of R. They are uniquely determined by R, and do not depend on the choice of 
fundamental cycles (i.e., they arc basis- independent). 

Substituting the Smith normal form (|3.11|) into (|3.10[) , we can write the condi- 
tions on O.2 as 

D aa <S> a = 27rm , 1 < a < g 2 , (3.12) 
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where m = P 1 n is integral and 3> = Qu>. $ a may be regarded as the flux of O2 
through a cycle C a in which the fundamental cycle appears with multiplicity 
given by Q a ^ 2 . 

From (|3.12[) . the allowed values of $ a depend on D aa , and in particular on 
whether D aa is one, greater than one, or zero. Let p denote the number of divisors 
of R equal to one, q = r —p the number of divisors greater than one, and s = <?2 — r 
the number of vanishing diagonal elements of D. Then (|3.12[) can be written as 

3>j = mod 27r, 1 < j < p, 

(f>p +fe = 2nm k /d k mod 27r, mfe = 0, . . . , dk — 1, 1 < k < q, 
$ r+/ = ai mod 2tt, < a t < 2tt, 1 < I < s, (3.13) 

where dk = D p+k ,p+k is shortened notation for the divisors of R greater than one. 
Once the <E> a 's are specified, u>, and hence f^, are determined by uj = 

Thus, topological gauge potentials are parameterized by s phases ai,...,a s 
taking values between and 2ir, and q phases 27rmi/c?i, . . . , 2-Km q /d q taking values 
constrained to be rational multiples of 2ir. We shall refer to the a/'s as free statis- 
tics phases and the 2irmk/ d^s as discrete statistics phases. Examples of both are 
given in Section 01 (More formally, the set of topological vector potentials modulo 
trivial gauge potentials, regarded as a group under matrix addition modulo 27r, is 
isomorphic to U(l) s x Z/di x •■■ x Z/d q , and is the character group of 7Ti(C2(r)).) 

(d) Aharonov-Bohm phases and two-body phases 

Among the free statistics phases, there may be some that correspond to one of 
the particles going around a cycle on G while the other remains fixed. Physically, 
such phases could be produced by solenoids threading the cycles of G with mag- 
netic flux (assuming the particles are charged). From the point of view of quantum 
statistics, we would like to distinguish between contributions to topological gauge 
potentials that may be attributed to individual particles interacting with an ex- 
ternal gauge potential, on the one hand, and contributions involving many-body 
effects on the other. 

Let fi be a gauge potential on a combinatorial graph G. We can construct a 
corresponding gauge potential on G2 , denoted fi^ 5 j following the prescription (|3.4[) 
for constructing a two-particle Hamiltonian from a single-particle Hamiltonian, i.e. 

^2-fj,Z}{fc.m} = fyk^lm + QjmO'kl + ^IkSjrn + ^hn^jk- (3-14) 

is antisymmetric (since il is antisymmetric), and hence constitutes a gauge 
potential on G2 ■ We shall call gauge potentials of the form (|3.14[) Aharonov-Bohm 
potentials. They are parameterized (up to a choice of gauge) by fluxes, or Aharonov- 
Bohm phases, 4>i B , ■ ■ ■ , (j>f B ' , through a set of fundamental cycles on G. 

Given a cycle C2 on G2, the flux Qf B (c2) can be evaluated as follows. Let 
p{c2) and q(c.2) denote the paths on G described by the individual particles as C2 
is traversed. If each particle returns to where it started, then p(c2) and 9(02) are 
themselves cycles, and we say that C2 is a direct cycle. In this case, Q^i ?) ^ s the 
sum of one-particle fluxes, i.e. 

ni B ( C2 ) = o(p( C2 )) + n( q (c 2 )). (3.15) 
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Note that (|3.15[) implies that Q 2 B is in fact a topological gauge potential; p(c £ .^,) 
and q(c t .<j,) are self-retracing cycles on the edges e and 4> through which the flux of fl 
obviously vanishes. If c 2 is not a direct cycle, then the particles exchange positions 
along C2, and we say that c 2 is an exchange cycle. In this case, the last vertex of 
p{c2) is the first vertex of 9(02), and vice versa, and 

n£ B (c 2 ) = n(p(c 2 )q(c 2 )). (3.16) 

We will say that two topological gauge potentials Q 2 and 2 are AB-equivalent 
if their difference, il 2 — Q,' 2l is an Aharonov-Bohm potential; cquivalently, f2 2 can be 
produced from f2 2 by adjusting the Aharonov-Bohm fluxes 4>f B - We would like to 
find parameters that determine topological gauge potentials up to AB equivalence. 
These parameters, together with the (f>^ B, s, then provide a complete parametriza- 
tion of topological gauge potentials (up to a choice of gauge) . 

We proceed by fixing a representative f2 2 from each family of AB-equivalent 
topological gauge potentials, and then finding parameters that characterize f2 2 . For 
the special case of circular graphs this is easily done. Circular graphs have their 
vertices arranged on a ring, and the only edges are between adjacent vertices. We 
may take f2 2 = 0, because every topological gauge potential f2 2 may be generated 
by an Aharonov-Bohm flux through the ring. Circular graphs arc discussed further 
in Section BE1 

Assuming that G is not a circular graph, we can choose a set of fundamental 
cycles Cj on G such that none of the Cj 's contains every vertex of G. Let kj denote 
a vertex of G not contained in Cj, and let c 2 -j denote the cycle on G 2 on which one 
particle traverses Cj while the other remains fixed at kj. We fix r2 2 ( U P to a choice 
of gauge) by requiring that 

fi*(c 2;i ) = mod 2tt (3.17) 

(a straightforward argument shows that every topological gauge potential is AB- 
cquivalent to a unique topological gauge potential satisfying (|3.17[0 . 

The / linear conditions (|3.17[) on f2 2 can be combined with the g 2 linear con- 
ditions (|3.6p which characterize general topological vector potentials, and the two 
sets of conditions written collectively as (cf. (|3.10[) ) 

R* -oj* = 2™*, (3.18) 

where R* is an integer matrix of dimension (/ + g 2 ) x f 2 , and n* is an integer 
vector of dimension (/ + 172)- The solution of (|3. 18[) proceeds as in Section l3fel 
Solutions are parameterized by s — / phases pi, ... , /3 S -/, which we call two-body 
phases, and q discrete statistics phases 27rmi/di, . . . , 2 r Km q /d q . It follows that a 
general topological gauge potential may be parameterized by / Aharonov-Bohm 
phases <pi B , . . . , 0^ B , which determine its fluxes through the C2 ; /s, in addition 
to the two-body phases /3; and discrete statistics phases 27rwfc/ 'dk- The lasso and 
bowtie graphs (Section 0EJ provide simple illustrations of the distinction between 
Aharonov-Bohm and two-body phases. 

(e) Bose and Fermi statistics 

Our treatment of quantum statistics on combinatorial graphs - referred to in 
what follows as the identified scheme - follows the approach of Leinaas and Myrhcim 



Article submitted to Royal Society 



14 



J.M. Harrison, J.P. Keating & J.M. Robbins 



(1977), treating the particles as classically indistinguishable. If instead one followed 
the conventional approach - referred to in what follows as the distinguished scheme, 
in which the particles are labeled from the start, one would find Bose and Fermi 
statistics. Below we establish the relationship between the two schemes. Note that, 
in the identified scheme, it does not make sense to speak of the exchange symmetry 
of a quantum state, as exchange is not defined - the wavefunction assumes a sin- 
gle value for each configuration of the indistinguishable particles. It turns out that 
Bose and Fermi statistics may be regarded as particular cases of the more general 
quantum statistics we have obtained here. In particular, Bose statistics corresponds 
to a trivial (e.g., vanishing) topological gauge potential, while Fermi statistics cor- 
responds to a topological gauge potential that assigns a phase of ir to exchange 
cycles and zero phase to direct cycles. 

Following the distinguished scheme one proceeds as follows. Let G denote a 
combinatorial graph with adjacency matrix A. We introduce the distinguished two- 
particle configuration space, denoted 02(G), or G 2 for short, that consists of ordered 
pairs j ^ I, of distinct vertices of G (note that particles are still not allowed 

to occupy the same vertex). We regard G2 as a graph with v(v — 1) vertices and 
adjacency matrix 

A2-jl,km = A ]k 5 lm + Al m Sjk- (3.19) 

Following the general quantization prescription given in Section QB we introduce 
the distinguished Hilbert space %2 = C v ^- V ~ 1 \ with basis vectors \jl), j ^ I, which 
describe states where particle 1 is at vertex j and particle 2 at vertex I. A quantum 
Hamiltonian on H2, denoted H2, is a v(v — l)-dimcnsional matrix with matrix 
elements 

H 2 . Akm = Ql\H 2 ftm). (3.20) 

satisfying (|2.4[) (so that short-time transitions are single-particle transitions). 

Exchange is defined on "H 2 by — > \lj). We decompose W.2 into two v(v — 
l)/2-dimcnsional subspaccs, H 2 , consisting of states that are even (a = 1) or odd 
(a = —1) under exchange. Indistinguishability is incorporated by requiring the 
Hamiltonian (and, indeed, any hcrmitian matrix representing an observable) to be 
invariant under exchange, i.e. H2-.ji.km = H2-,ij, m k- Then H2 preserves exchange 
symmetry, i.e. it leaves the subspaces W invariant. 

For example, given a one-particle Hamiltonian H , we may construct the two- 
particle Hamiltonian 

H 2 -jLkm = HjkSlm + Hi m Sjk- (3-21) 

On the antisymmetric subspace T~L 2 , the energy levels of H2 are just sums of distinct 
energy levels of H , and the two-particle eigenstates are antisymmetric products of 
one-particle eigenstates - this corresponds to noninteracting fermions. The situation 
is different on the symmetric subspace H 2 ■ In general, the spectrum of H2 on H 2 
is not simply related to the one-particle spectrum (although in some special cases 
it is). This is because the diagonal states \jj) are excluded; bosons are prevented 
from occupying the same vertex, which, unlike fermions, they would otherwise do. 

The equivalence between the identified and distinguished schemes may be es- 
tablished by means of a unitary map, U a , from the identified Hilbert space, T-L 2 , to 
the symmetric or antisymmetric subspace, 7-L 2 j °f the distinguished Hilbert space. 
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U a is given by 

U°\jl) = -^=(\Tl) + o-\Tj)), j<l. (3.22) 

Given an exchange- invariant Hamiltonian H 2 on , we can construct a unitarily 
equivalent Hamiltonian H 2 on ^2 by taking = U a H 2 U a ^ . The matrix elements 
of iff are given explicitly by 

H 2;jl,km = H 2 -jl,km + <jH 2 -lj,k m , j <l,k < 171. (3.23) 

Thus, in the identified scheme, the statistics sign a appears in the Hamiltonian H% , 
while the Hilbert space, ri 2 , is independent of a. In contrast, in the distinguished 
scheme, a determines the Hilbert space ri 2 > while the Hamiltonian H 2 is indepen- 
dent of a. In particular, if H 2 is given by the sum of one-particle Hamiltonians as 
in (|3~2"Tj) . then (|3~2^|) yields 

(jl\H 2 \km) = 5 jk H lm + aS-jmHik + aSi k H jm + Si m H jk , j <l,k < m. (3.24) 

This is in accord with (|3.4p . and establishes that a corresponds to Bose or Fermi 
statistics. 

From (|3.23p and the constraints (|2.4[) . we can express the relation between the 
Bose and Fermi Hamiltonians H 2 an d ^2 as 

H 2 - n , km = e ia ^H+ ml (3.25) 



where Cl 2 is given by 



7r, j = m,l ~ k, 

^2;.jl,km = { j ~ m ! ' = ( 3 ' 26 ) 

0, otherwise. 

Thus, H+ and # 2 ~ are related by a gauge potential, as fl 2 is real antisymmetric and 
its (j7,fcm)th element vanishes unless {j,l} ~ {fc,m}. In fact, fif" is a topological 
gauge potential. This can be verified by checking the conditions (|3.6p explicitly. 
Alternatively, we may argue as follows: given a path p 2 on G 2 of length 7-, let 
p(P2) = {xi, ■ ■ ■ ,x r ) and q 2 (p) = (y, ■ ■ ■ , y r ) denote the single-particle paths along 
p 2 . Nonzero contributions of ±7r to fl 2 (p 2 ) come from those edges of p 2 where 
the ordering of the single-particle positions switches over, either from Xj < yj to 
Xj+i > Uj+i or from Xj > yj to Xj+\ < yj+i- However, on a direct cycle c 2 , and, 
in particular, on c e ^, the particles separately return to where they started, so the 
number of switchovers, and hence the number of ±7r contributions, must be even - 
therefore Q, 2 {c 2 ) is a multiple of 2tt. Indeed, by this reasoning we see that 



1, if c 2 is a direct cycle, ^ ^ 

i—l, if c 2 is an exchange cycle. 

4. Examples 

We investigate statistics phases for a number of graphs, starting with the simplest 
examples. 
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(a) Linear graphs 

The linear graph Ljv consists of N vertices on a line with adjacent vertices 
connected by edges, so that Aj^ = <Ju_fcM. From the point of view of quantum 
statistics, linear graphs are trivial; it turns out that there are no nontrivial topolog- 
ical gauge potentials on two-particle linear graphs. However, linear graphs provide 
simple examples where the free-particle energy levels and eigenstates can be calcu- 
lated explicitly, and they serve to illustrate some points of the preceding discussion, 
including the equivalence of Bose and Fermi statistics in cases where topological 
phases are absent, as well as the effect of non-topological gauge potentials. 

Let H be the one- particle kinetic energy on Lm (cf. Section It is straightfor- 
ward to show that the energy levels of H are given by E a = 4 sin (ixa/2N), < a < 
N, while the eigenstates are given, up to normalization, by (j\ip a ) = cos(naj/N). 
We take the two-particle Hamiltonian, H^, to be the sum of one-particle Hamil- 
tonians as given by (|3.4|) . As discussed in Section 0j!J a = 1 corresponds to Bosc 
statistics and a = — 1 to Fermi statistics. It is straightforward to show that the 
energy levels of are independent of a, and are given by sums E a + Ef, with a 
and b distinct. The corresponding eigenstates are, for a = — 1, antisymmetric prod- 
ucts of distinct eigenstates \ip a ) and \ipb), and, for a = +1, symmetric products of 
these eigenstates with diagonal components removed. The interaction between the 
(hard-core) bosons is reflected in the fact that two-particle states with a = b are 
absent from the spectrum. 

The fact that the Bose and Fermi spectra coincide for linear graphs can be un- 
derstood from our topological treatment of quantum statistics. The only nontrivial 
cycles on the two-particle configuration space correspond to each particle moving 
back and forth in alternation along separated segments of the linear graph. Such 
cycles are metrically contractible, so that there are no nontrivial topological gauge 
potentials. It follows that Fermi gauge potential f2f (cf. p.26p ) is just a gauge 
transformation. 

There is an alternative, single-particle interpretation of H% , namely as a discrete 
approximation to the quantum Hamiltonian for a free particle in a two-dimensional 
right triangular domain with sides 1, 1 and v2 (i.e., half the unit square below the 
diagonal), and with h = 1/N. Neumann boundary conditions apply on the right 
sides of the triangle, while a determines the boundary condition on the diagonal 
(Neumann for a = 1, Dirichlct for a = — 1). 

Fixing the Hamiltonian to be for definiteness, it is instructive to observe the 
consequences of introducing a gauge potential, denoted (l 2 , that is not topological. 
Let us take to have flux 2irp/t through every cycle c £i where e is an edge between 
vertices in the range between r and r+t, and cf> an edge between vertices in the range 
between s and s+t. For cycles outside this range, we take f22(c e .0) = 0. In terms 
of the single-particle interpretation, 0,2 corresponds to a uniform magnetic field of 
strength 2-Kp/{t/N) (in units where q/c, q being the charge of the particle, is equal 
to one) through the square of area t 2 /N 2 with diagonal corners (r/N,s/N) and 
((r + t)/N, (s + t)/N). Numerical calculations show that the shifts in energy levels 
produced by ^2 scale with 1/N, in accord with the fact that a magnetic field does 
not change the mean density of states. However, with VL2 present, while typical 



eigenstates are delocalized (cf. Figure 2(a) I, one finds some eigenstates that are 



strongly localized in the flux square 10 < j < 15 and 25 < / < 30 (cf. Figure 2(b) ) 
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In the single-particle interpretation these are Landau-like levels, and are indicative 
of Lorentz forces in the classical dynamics. 



Typical eigenstate 



Landau-like level 
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(b) 



Figure 2. Consecutive eigenstates for two particles on a linear graph with N = 40 vertices 
with gauge potential (a) Typical eigenstate (b) Eigenstate localized in the flux square, an 
analogue of a Landau-level wavefunction 



(6) Circular graphs 

The circular graph, Cat, is obtained by connecting the first and last vertices of 
the linear graph Ljy. For a circular graph, we naturally expect, and indeed find, 
anyon statistics, confirming that our model provides a reasonable description of 
quantum mechanics on a loop. For simplicity we consider a loop with three vertices, 
shown in Figure |3] (a) (the conclusions are similar for N > 3). Writing down the 
three allowed two- particle configurations in Figure [3] (b), it is apparent that the 
two-particle graph G2 is also a single loop. A single traversal of this loop is an 
exchange cycle, and the associated flux <j) corresponds to anyon statistics. The flux 
4> can be generated by an Aharonov-Bohm potential, corresponding to a solenoid 
threading the one-particle graph C3. 



a) 




1 3 (1,2) (2,3) 

Figure 3. (a) C3 (b) Two-particle graph for C3 



(c) Star graphs 

The star graph, S e , shown in Figure 0] (c), consists of e vertices each connected 
to a central vertex, and so has e edges and e + 1 vertices. We consider first the 
e = 3 star graph, or y-graph, for which the two-particle graph is easily displayed 
(Figures H] (a) and (b)). The two-particle graph consists of a single cycle which 
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exchanges the particles through the arms of the 'Y'. A flux through this cycle 
produces anyon statistics. For e > 3, the two-particle star graph consists of i>2 = 
(e + l)e/2 vertices and &2 = e(e — 1) edges. The number of independent cycles, fa, 
is given by (e — l)(e — 2)/2. There are no Aharonov-Bohm phases (there are no 
nontrivial cycles on S e ) and no constraints on topological gauge potentials (there 
are no disjoint edges on S e ). Therefore, topological gauge potentials on S e are 
parameterized by gi = (e — l)(e — 2)/2 two-body statistics phases. The number 
of statistics phases can also be obtained from the following simple argument: each 
phase corresponds to the choice of a pair of edges along which to exchange the 
particles, given that the particles start on the vertices of some given edge. g% is 
therefore the number of pairs of e — 1 objects. 



2 (1,2) (2,4) 3 




a) b) c) 

Figure 4. (a) Y-graph (b) Two-particle Y-graph (c) Star graph with e edges 



(d) Lasso and bowtie 

The lasso graph consists of a three- vertex loop with a single external lead - see 
Figure [5] (a) . It provides a simple example of quantum statistics which combines 
aspects of circular graphs and star graphs discussed above. The two-particle lasso, 
G2, is shown Figure[5](b) along with a spanning tree, T 2 ( m bold). We fix the gauge 
so that the edges of T2 are assigned zero phase. The three edges of G2 not in T2 
determine the fundamental cycles. The central square corresponds to the metrically 
contractible cycle c £i¥ , on which the particles move in alternation along the edges e 
and ip of the lasso. As the flux through this cycle must vanish, the edge ((1, 3), (2, 3)) 
is assigned zero phase as well. The left triangle of G2 corresponds to the cycle in 
which the one of the particles goes around the loop of the lasso while the other 
remains on the external lead, and may be assigned an Aharonov-Bohm phase (f> AB . 
The right triangle of G2 corresponds to an exchange cycle in which the particles 
move around the loop of the lasso. The associated phase, denoted a, is a two-body 
phase. The length-six cycle along the perimeter of G2 coincides with the exchange 
cycle on the Y graph, and has phase a + 4>ab- These results coincide with those 
of Balachandran & Ercolessi (1992), who considered the two-particle metric lasso 
graph. 

Another example considered by Balachandran & Ercolessi (1992) is the bowtie, 
or figure-of-eight, which consists of two three- vertex loops sharing a common vertex 
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- see Figure[3](c). Calculations show that the two-particle bowtie has two Aharonov- 
Bohm phases (corresponding to the two loops) and two two-body phases. 



(1,4) (2,4) 




Figure 5. (a) Lasso (b) Two-particle lasso (c) Bowtie 



(e) Nonplanar graphs: K3 3 and the molecule 

K 5 is the fully connected graph with five vertices, and K 3:3 is the fully connected 
bi-partite graph with two sets of three vertices - see Figure[6](a) and (b). A theorem 
of Kuratowski (1930) states that every non-planar graph (i.e., a graph that cannot 
be drawn in the plane without crossings) contains K 5 or AT 3:3 as a subgraph (pos- 
sibly after contracting some edges to points) . From the point of view of quantum 
statistics, K§ and A33 are interesting because they are the smallest graphs that 
exhibit discrete statistics phases. Calculations following the procedure of Section I5fel 
(details are omitted) show that K$ has six Aharonov-Bohm phases and ^3.3 has 
four Aharonov-Bohm phases (corresponding to their respective number of funda- 
mental cycles). In addition, both have a single discrete statistics phase that can be 
either or ir (mod 2ir). Cycles whose fluxes are given by this discrete phase are 
necessarily exchange cycles. An example for K$ consists of the cycle in which one 
particle goes around a three-vertex loop (e.g., (1,2,3,1)) while the other remains 
fixed, followed by an exchange around the same three- vertex loop in the opposite 
direction (e.g., ((1, 2), (2, 3), (1, 3), (1, 2))). An example for A' 3j3 consists of the cy- 
cle in which one particle goes around a four-vertex bowtie loop (e.g., (1,4, 2, 5, 1)) 
while the other remains fixed, followed by an exchange around the same four-vertex 
loop in the opposite direction (e.g., ((1,4), (4, 5), (1, 5), (1, 2), (1,4))). 

The K 5 molecule is the graph consisting of two K5S joined by a single edge, 
as in Figure [5] (c) . Calculations (again, details are omitted) show that two discrete 
statistics phases appear, both either or n, along with 12 Aharonov-Bohm phases 
and 6 two-body phases. 

5. Discussion 

The abelian statistics of two indistinguishable quantum particles on a combina- 
torial graph are characterized by a set of continuous and discrete- valued phases. 
The continuous phases may be separated into some that are produced by external 
Aharonov-Bohm fluxes and the rest describing two-body statistical interactions. 
The appearance of discrete phases may be related to whether the graph is planar 
or not, a connection that merits further investigation. While we have concentrated 
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on the case of two particles, the abelian statistics for more than two particles may 
be characterized and calculated using the results presented here. 

Nonabelian statistics requires new considerations. A full description is needed 
of the fundamental group of the n-particle configuration space, or, equivalcntly, 
the braid groups of the graph (not simply their abelianized versions). Recently, 
Farley & Sabalka (2005) have developed methods based on discrete Morse theory for 
obtaining efficient presentations of graph braid groups. We will discuss applications 
of these techniques to nonabelian graph statistics in a forthcoming publication 
(see (Kitaev 2008) for an exact solution of a spin-lattice model where nonabelian 
statistics emerge). 

From the point of view of physics, one of the principal attractions of quantum 
graphs is that they provide mathematically tractable models of complex physical 
systems. As applications have so far concentrated on independent-particle models, 
the scope for manifestations of quantum statistics on graphs is great. We suggest 
a few possibilities here. Graph statistics may play a role in many-electron net- 
work models of molecules, in analogy with the emergence of the Berry phase - the 
molecular Aharonov-Bohm effect (Mead and Truhlar 1979) - in molecular spectra. 
Topological signatures in single-particle transport on networks (Avron 1995), which 
provide models and variants of the quantum Hall effect, may have many-particle 
generalizations in which statistics plays a role. Many-particle graphs may provide 
new models for anyon superconductivity (Wilczek 1990). An intriguing potential ap- 
plication of nonabelian graph statistics is to topological quantum computing (Nayak 
et al. 2008). There one looks for systems with a degenerate ground state spanned 
by distinct quasiparticle configurations, in which the only (easily) realizable evolu- 
tions are, up to a phase, a discrete set of unitary transformations generated by the 
(adiabatic) exchange of quasiparticlcs. By introducing spin (Harrison 2008), quan- 
tum graphs might also provide models in which to investigate the role of quantum 
statistics in the quantum spin Hall effect and topological insulators (Hasan & Kane 
2010, Qi & Zheng 2010). 

Applications will depend on nontrivial graph statistics emerging in a particular 
model. A standard approach would be to look for novel many-particle ground states 
on graphs and to study their excitations. This is work for the future. However, even 
without a specific mechanism, we believe it is worthwhile to pursue a general inves- 
tigation of quantum statistics on networks and its consequences. Quantum physics 
has provided many examples where topology underlies new and unexpected phe- 
nomena. Nature seems to exploit the opportunities available to it, and discoveries 
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may follow from knowing where to look and what to look for. 
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